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Statistical  OKuacterization  of  Rough  Terrain 


1.  INTRODUCTION 

In  this  report,  a  general  technique  for  characterizing  data  will  be  developed. 

To  illustrate  the  procedure,  a  particular  application  to  terrain  features  will  be  used 
as  an  example. 

The  general  problem  we  are  formulating  is  the  following.  We  consider  a  set 
of  data  and  whatever  external  requirements,  results,  experience,  and  limitations 
are  available,  and  propose  statistical  distributions  which  might  characterize  the 
data.  Next,  we  mathematically  derive  a  multivariate  joint  probability  dersity  func¬ 
tion,  PDF,  for  each  case.  W*»  then  use  statistical  estimation  theory  to  evaluate 
the  parameters  of  the  PDF's,  and  finally,  we  apply  an  hypothesis  test  to  decide 
which  is  the  more  appropriate  PDF  for  the  data  set  or  subsets.  This  general  pro¬ 
cedure  involves  a  number  of  topics  and  some  further  discussion  may  clarify  these 
aspects. 

Consider  some  arbitrary,  randomly  distributed  physical  quantities: 

Zjt  Zg,  ...,  Z^— either  discrete  or  continuous.  One  set  of  measurements  of  these 

physical  quantities,  consisting  of  one  value  of  each  of  N  variables:  z,,  z„,  ....  z„ 

12  N 

is  considered  as  a  single  representation  of  a  random  process  governed  by  some 
multivariate  PDF.  The  PDF  can  have  different  forms,  Gaussian,  exponential. 
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uniform,  binomial,  and  so  on.  Properties  of  the  multivariate  PDF  are  discussed 

1  2 

in  numerous  texts,  such  as  Papoulis  or  Mood  and  Graybill.  In  general,  the  PDF 

is  a  function  of  several  parameters  such  as  the  mean  values  of  each  of  the  variables 
2 

u.,  the  variances  o. ,  and  the  covariances  y..»  where  i  =  1,2,  ....  N,  and 
'l  i  ij 

j  =  1,  2,  ....  N. 

When  starting  from  observed  values,  the  situation  in  its  most  general  state 
would  be  that  neither  the  form  of  the  PDF  is  known,  nor  the  parameters  which  enter 
into  the  expression  for  the  PDF.  In  order  to  obtain  some  information  on  these 
parameters,  estimation  theory  is  used.  The  next  question  is  that  of  the  form  of 
the  PDF  in  which  the  parameters  are  incorporated.  This  is  not  known  a  priori  and 
must  be  determined.  An  hypothesis  testing  procedure  is  developed  to  accomplish  this. 

The  particular  hypothesis  testing  procedure  selected  in  the  present  case  allows 
only  a  binary  decision  process  to  be  considered.  Hence,  the  discrimination  is  re¬ 
stricted  to  two  forms  of  the  PDF.  The  test  is  based  upon  the  maximum  a  posteriori 
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probability  criterion.  This  is  equivalent  to  the  minimum  error  probability  criter¬ 
ion.  The  explicit  mathematical  expressions  involved  in  the  hypothesis  test  are 
worked  out  for  the  case  of  choosing  whether  an  N  -  variate  Gaussian  or  exponential 
PDF  would  be  more  likely  to  produce  the  occurrence  of  the  data  set  [z^].  However, 
the  procedure  for  hypothesis  testing  may  be  applied  quite  generally.  Indeed,  w-here 
the  number  of  cases  allow,  a  better  test  would  be  to  see  which  density  would  be 
more  likely  to  result  in  several  realizations  rather  than  the  single  one  used  here. 

The  specific  problem  used  as  an  example  in  this  report  is  that  of  characterizing 
a  large  terrain  region  that  is  considered  to  be  made  up  of  smaller  subareas 
(~4  km  ).  The  main  feature  of  interest  is  the  distribution  of  heights  within  these 
subregions.  The  nature  of  this  specific  problem  places  some  constraints  on  the 
form  of  the  particular  statistical  analysis  carried  out  for  the  report. 

The  goal  of  the  particular  case  is  the  development  of  mathematical  descriptions 
of  each  of  the  subareas  for  use  in  calculation  of  the  scattering  of  electromagnetic 
waves  from  the  uneven  terrain  surface.  Each  region  is  characterized  by  a  geologic 
code  and  several  statistical  parameters.  In  particular,  we  are  concerned  with  being 
able  to  associate  a  PDF  with  the  range  of  heights  in  the  subregions  and  to  deter¬ 
mine  parameters  that  make  the  general  PDF  explicit. 


1.  Papoulis,  A.  (1965)  Probability,  Random  Variables  and  Stochastic  Processes, 

McGraw-Hill. 

• 

2.  Mood,  A.M. ,  and  Graybill,  F.A.  (1963)  Introduction  to  the  Theory  of  Statistics, 

McG raw  Hill. 

3.  Jenkins,  G.  M. ,  and  Watts,  D.G.  (1968)  Spectral  Analysis  and  Its  Applications, 

Holden-Day. 

4.  Whalen,  A.  D.  (1971)  Detection  of  Signals  in  Noise,  Academic  Press. 
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The  data  set  [2.]  which  is  analyzed  in  detail  in  this  report  is  such  that  each 
value  of  the  random  variable  corresponds  to  a  point  in  the  x  -  y  plane,  so  that 
z.  =  z.  ,xk,  y^),  where  i  =  1,  2,  3,  . . .  N,  and  N  is  the  total  number  of  grid  points 
in  the  x  -  y  plane.  Here,  denotes  the  k*h  equally  spaced  x-value  along  the 
x-axis  and  y^  denotes  the  £lh  y-value  along  the  y-axis,  where  k  =  1,  2,  .... 

and  £  =  1,  2 .  ^W.  Thus,  the  N  points  are  dis  ributed  in  the  x  -  y  plane  so 

as  to  form  a  rectangular  grid.  This  restricts  the  analysis  to  the  class  of  data  sets 
where  the  covariance  matrix 


R  „  =  <(z  -  z)  (z  -  z)  > 

mn  m  n 


can  be  assumed  to  have  the  form: 


R  =  a2  exp(-r2  /T2) 
mn  mn 


(1) 


(2) 


where 

2  _  _ 

a  *  variance  <(z  -z)(z  -  z }  > 

P  P 

z  =  mean  value 


<  >  =  denotes  expectation  value 

T  =  correlation  length 


and 


2 

Tmn 


(x 


m 


2  2 
x  )  4-  (y  -  y  )  that  is, 
n  m  n 


=  square  of  distance  between  grid  point 

(x  ,  y  )  and  point  <x  .  y  )  . 
m  m  n  n 


The  motivation  for  assuming  a  covariance  matrix  in  the  form  given  by  Eq.  (2) 
arises  because  of  the  fact  that  it  leads  to  a  tractable  mathematical  expression  for 
the  incoherent  power  that  is  scattered  when  an  electromagnetic  wave  is  reflected 


from  a  rough  surface.  '  .Then,  the  data  set  [Zj]  corresponds  to  the  heights  at 
particular  points  on  a  rough  surface. 

Also,  in  the  hypothesis  testing  procedure  the  binary  decision  process  is  worked 

out  in  detail  here  for  the  case  where  the  PDF  is  either  an  N  -  variable  Gaussian  or 

exponential.  This  specialization  is  also  motivated  by  the  theory  of  electromagnetic 
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wave  scattering  from  rough  surfaces.  * 

The  performance  of  a  radar  system  in  detecting  and  tracking  targets  depends 

upon  the  electromagnetic  scattering  characteristics  of  the  terrain  surrounding  the 

radar.  Rough  terrain  will  contribute  to  the  radar  clutter  return  and  to  the  multi- 
5  6 

path  return  *  which  involves  wave  scattering  from  both  the  target  and  terrain 
surrounding  the  radar.  These  two  aspects  can  be  described  by  the  theory  of  scat¬ 
tering  from  rough  surfaces  if  the  properties  of  the  terrain  surrounding  the  radar 
are  known.  The  properties  of  the  terrain  pertinent  to  rough  surface  scattering 
include  the  form  of  the  PDF  for  the  surface  height  distribution,  the  correlation 
length  when  Rmn  is  assumed  to  have  the  form  given  by  Eq.  (2),  the  mean  surface 
height,  the  variance  in  surface  height,  and  the  complex  dielectric  constant 
characterizing  the  surface. 

In  this  introduction,  we  have  proposed  a  general  technique  for  statistical 
characterization  and  discussed  how  it  has  been  applied  to  a  specific  case.  The 
particular  external  considerations  used  in  defining  appropriate  forms  and  constraints 
for  the  formulation  have  been  included.  In  the  main  body  of  this  report  we  will 
develop  the  successive  steps  used  in  the  analysis  as  exemplified  by  the  terrain 
height  data  characterization.  The  first  aspect  is  the  parameter  estimation  required 
for  the  eventual  selection  of  an  appropriate  PDF. 


2.  PARAMETER  ESTIMATION 

The  starting  point  in  these  procedures  is  a  data  set  or  sets  w!u  :h  are  assumed 
to  represent  samples  characterized  by  appropriate  PDF’s.  For  t’  e  example  with 
which  we  are  concerned  our  potential  PDF's  are  Gaussian  or  e  .ponential.  The  first 
step  in  discriminating  between  these  is  to  construct  estimators  or  the  parameters 
required  for  the  two  different  densities.  The  estimators  can  be  considered  as  func¬ 
tions  defined  on  the  N-dimensional  sample  space.  The  particular  estimate  obtained 
from  the  given  data  is  then  regarded  as  a  realization  of  the  random  variables 


Beckmann,  P. ,  and  Spizzichino,  A.  (1963)  The  Scattering  of  Electromagnetic 
Waves  from  Rough  Surfaces.  Macmillan  Co. 

Ruck,  G.T. ,  Barrick,  D.  E. ,  Stuart,  W.  D. ,  and  Krichbaum,  C.  K.  (1970) 
Radar  Cross  Section  Handbook,  Vol.  2,  Plenum  Press. 
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representing  the  estimator.  The  parameters  which  will  be  discussed  are  the 
means  p .,  the  variances  cr.  ,  and  the  covariance  function  for  the  two  possible 
PDF's, 

For  a  completely  rigorous  approach  to  estimation,  we  could  use  a  technique 
such  as  maximum  likelihood  estimation.  That  estimate  of  the  parameter  maxi¬ 
mizes  the  probability  of  having  obtained  the  observed  values.  However,  when  the 

data  represents  a  multivariate  correlated  distribution,  the  procedures  for  maxi- 

7  .  . 

mizing  the  likelihood  functions  become  complex  and  require  multiple  observations 
of  the  N -variate  samples.  As  an  alternative,  we  would  consider  the  accuracy  of 
the  estimator  by  minimizing  the  bias  and/or  variance  of  the  estimator,  which  are 
defined  in  terms  of  the  lower  order  moments  of  the  estimator  PDF. 

The  bias  of  an  estimator  is  defined  as  the  difference  between  the  true  value  of 
the  parameter  and  the  expected  value  of  the  estimator.  When  the  bias  is  zero,  the 
PDF  of  the  estimator  is  centered  exactly  at  the  true  value  of  the  parameter,  and 
the  estimator  is  termed  unbiased.  This  implies  that  as  the  number  of  observations 
increases,  the  estimate  of  a  parameter  approaches  the  true  value  of  the  parameter. 
Also,  if  the  variance  of  an  estimator  is  small,  the  PDF  of  the  estimator  has  a 
narrower  width  centered  around  its  mean  value.  If  the  bias  and  variance  of  an 
estimator  both  tend  to  zero  as  the  number  of  observations  increases,  then  the  esti¬ 
mator  is  termed  consistent.  However,  not  all  estimators  are  consistent,  since 
reducing  the  variance  may  sometimes  Increase  the  bias.  A  compromise  between 
minimum  variance  and  minimum  bias  is  often  achieved  by  minimizing  the  mean 
square  error  of  the  estimator.  The  problem  with  this,  however,  is  that  assignment 
of  values  to  the  bias  and  variance  of  the  estimator  depends  on  a  knowledge  of  the 
PDF  of  the  estimator  which  in  turn  is  related  to  the  PDF  of  the  variates.  In  general, 
unless  the  variates  are  independent  with  a  Gaussian  PDF,  it  is  not  easy  to  establish 
the  appropriate  PDF's  for  the  parameter  estimators  and  hence  the  bias  and  variance 
cannot  readily  be  evaluated. 

For  a  specific  problem,  that  is,  a  single  terrain  region  which  we  describe  by 
a  set  of  N  height  values,  it  would  be  possible  to  carry  out  p-repeated  measurements 

[z.Jj, _ _  [z.J  p  of  the  N -variates  at  appropriate  positions  and  then  do  a  maximrr., 

likelihood  estimate  (mle)  of  the  parameters  along  the  lines  suggested  by  Morrison* 
for  a  multinormal  case.  If  his  analysis  is  extended  to  our  multivariate  exponential 
case,  though,  we  see  that  the  mle  for  the  mean  and  the  covariance  (and  variance) 
estimators  are  not  as  straightforward.  Both  forms  would  have  to  be  calculated  for 
consistency.  In  our  actual  case,  there  are  a  multitude  of  regions  each  wit*  differ¬ 
ent  density  functions;  then  the  full  mle  analysis  would  have  to  be  done  for  each  case 
and  again  would  become  computationally  unappealing. 

7,  Morrison,  D.  F.  (1976)  Multivariate  Statistical  Methods,  McGraw-Hill. 
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This  discussion  has  pointed  out  why  we  did  not  conduct  an  assessment  of  the 
accuracy  of  the  parameter  estimators  in  any  strict  sense.  This  may  not  be  the 
case  for  other  situations  and  based  on  the  particular  circumstances,  some  assess¬ 
ment  of  the  confidence  limits  and  error  bounds  may  cause  the  user  to  select  some 
more  appropriate  forms  for  his  estimators. 


For  our  case  of  terrain  heights,  the  number  of  subareas  to  be  evaluated  and 
the  difficulty  of  obtaining  successive  measurements  of  our  N-variates  that  would 


essentially  be  independent  observations  (given  the  digitized  terrain  maps  available 
for  analysis)  led  us  to  the  use  of  a  single  N-variate  sample.  We  thus  had  to  make 
a  number  of  assumptions.  We  have  assumed  that  the  mean  height  for  each  variate 
is  the  same  and  the  variances  of  each- of  the  heights  about  their  means  also  are 
equal.  With  these  assumptions  it  becomes  reasonable  to  use  the  sample  mean  as 
the  estimator  for  the  mean  height: 


£  Z(x.,  y.) 
j=l  1  J 


and  as  the  estimator  for  the  variance,  the  sample  variance: 

oz  =  u m  £  £  (z..  -z)2 

i=i  j=l 


where  the  mean  is  given  by  its  estimated  value. 

In  constructing  an  estimator  for  the  covariance  matrix  parameter  ~‘7Z> 
process  is  slightly  more  complicated.  We  assume  that  the  height  distribution  is 
stationary  over  the  individual  subregions,  although  the  relations  for  distinct  areas 
may  differ.  This  is  a  spatial  equivalent  to  a  time  series  that  is  stationary  for  some 
interval  although  not  over  a  very  long  duration.  This  again  requires  the  means  and 
variances  of  the  PDF's  to  be  independent  of  location  in  the  subregion  and  the 
covariances  to  be  dependent  only  on  a  discrete  separation  factor  for  all  spatial  posi¬ 
tions.  This  concept  of  position  being  introduced  into  the  covariance  relations  has 
been  discussed  in  the  introduction  but  there  are  a  number  of  aspects  to  this  depend¬ 
ence  which  affect  toe  form  of  the  estimator,  and  hence  should  be  given  further 
attention. 

The  objective  of  our  covariance  analysis  is  to  generate  a  functional  form  for 

the  covariance  between  variates  at  two  points  of  the  grid  structure  that  is  based 

solely  on  their  separation  (that  is,  the  process  is  stationary).  In  particular  we 

define  a  correlation  length  T,  where  T  represents  the  distance  at  which  a  normalized 

covariance  function,  t  has  decreased  to  the  value  e  1  where  e  *  .  Based 

zz  zz  a 
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on  the  grid,  c  takes  on  discrete  values  and  c  (0)  =  1.  For  the  calculation  of 
°  zz  zz 

T,  we  assume  a  parabolic  fall-off  with  distance  and  carry  out  a  least  squares  fit 
to  the  function.  Separate  x  and  y  direction  calculations  were  made  and  averaged 
for  uniformity.  In  addition,  the  fall-off  is  rapid  enough  so  the  least  squares  fits 
were  based  on  averaging  the  first  x  and  y  separations  respectively. 

One  implication  of  this  procedure  is  that  direct  estimation  of  the  entire  covari¬ 
ance  matrix  was  not  carried  out.  As  indicated,  the  short  -’ange  covariances  were 
used  to  determine  a  \alue  for  T  and  subsequently  the  complete  covariance  matrix  was 

-T  2  /T2 

*  A  2  1  mn  x 

estimated  by  appropriate  functional  dependence  y  ( m,  n)  =  <r  e  .  If  we 

had  not  been  interested  in  writing  the  covariance  matrix  in  functional  form  we  could 

have  chosen  the  alternative  approach  of  specifying  y  '  (r )  where  as  before 
2  2  2 

i  =  (x  -  x  )  +  (y  -  y  )  -for  every  pair  of  variates  in  the  terrain  region, 
mn  '  m  n  vm  ■'n  J  r  " 

The  actual  forms  used  are: 

Vn  rVN  V^-k  1  *0 

e‘<tt=,1,N52,[§  5, 


tVn  ■Jii-k  .[Vn  Vn-k 

£  E  hjhi^-^E  E  hr 


i=l  j=l 


i=l  j=k+l 


Under  the  assumptions  that  have  been  made,  including  stationarity  over  the  region, 
the  above  forms  for  the  estimators  are  similar  to  the  form  used  in  Jenkins  and 
Watts.  3 

The  least  squares  results  for  these  values  give  us  correlation  lengths  for  the 
x  and  y  directions  and  their  average  is  used  as  the  value  of  T  that  applies  to  the 
entire  region.  This  then  generates  the  complete  covariance  matrix. 

The  whole  question  of  the  estimation  theory  used  in  the  particular  case  dis¬ 
cussed  in  the  report  is  complex.  Care  should  be  exercised  by  other  users  of  this 
approach  to  apply  estimators  which  are  appropriate  to  their  knowledge  of  the  data 
they  are  analyzing.  In  the  present  case  we  are  in  the  process  of  considering  the  use 
of  multiple  observations  to  establish  better  estimators  for  the  parameters  of  the 
probability  density.  A  further  question  of  interest  is  the  relation  of  the  estimators 
to  the  outcome  of  the  hypothesis  testing  procedure  in  which  they  are  used.  This 
aspect  is  also  being  addressed. 

In  addition  to  the  estimators,  the  hypothesis  testing  also  requires  the  forms  of 
the  PDF's  being  evaluated.  The  derivation  of  the  two  PDF's  used  in  the  present 
example  constitutes  the  next  topic. 


11 


3.  MULTIVARIATE  EXPONENTIAL  PROBABILITY  DENSITY 


In  our  analysis  of  the  terrain  characteristics,  we  want  to  establish  wMther  the 
heights  within  the  selected  regions  are  more  closely  distributed  according  to 
Gaussian  or  exponential  probabilities.  As  the  first  step  in  this  determination  we 
selected  estimators  for  the  mean  and  variance  of  the  distribution  to  be  verified. 
Here  we  establish  the  general  N  -  variate  forms  for  the  two  PDF's  in  which  we  are 
interested.  We  are  assuming  that  the  N  -  variates  are  jointly  distributed  and  all 
have  the  same  mean  and  variance. 

2 

The  form  for  the  Gaussian  multivariate  PDF  is  well  known.  Let 
2.'  =  (z  z2  •  •  •  zjJ)  be  an  N  -  dimensional  random  variable  in  vector  form.  Then 

Pg<4  4  -  Z-)  .  (|S|'/W-)-V(1-<5'  •  t,Tf1  (£'  -  (61 


where  R  is  a  positive  definite  symmetric  matrix  whose  elements  are  constants  and 

fi  is  a  vector  fi  =  fi g*  Mj^)  whose  elements  are  constants. 

For  our  case  fi.  =  fi,  Vi  and  the  matrix  R  is  the  covariance  matrix  for  the 
1  2  ~  i 

variates  where  R;,  =  O;.,  O;;  =  ct  ,  Viandp..  =-%*-  .  Thus 

1J  11  ”  (T 4 


R  =  or 


PN1  PN2 


Note  that  the  quadratic  form  for  this  case  is  given  by 


Q2  =  (z'  -  fi)1  R"1  (z'  -  fi) 


To  complete  the  analysis  it  is  still  necessary  to  derive  an  equivalent  form  for 
a  multivariate  exponential  PDF.  In  that  case  the  appropriate  coefficients  will  be 
required.  The  general  form  would  be 


PE  (zl’  z2’  ZN}  =Ci® 


-C2[(z*  -  fi)T  R'V'  -  /£)] 
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where  the  quadratic  form  is  defined  as  in  the  case  for  the  Gaussian  PDF.  The 
question  to  be  resolved  before  this  can  be  used  is  what  values  must  be  assigned  to 
Cj  and  C2- 

The  values  can  be  determined  by  evaluating  the  zeroth  and  second  moment 
integrals.  For  a  probability  density,  the  zeroth  case  must  equal  unity  and  the 
second  moment  equals  the  variance.  The  two  integrals  turn  out  to  be  similar  in 
man\  aspects  of  their  evaluation.  We  first  discuss  the  zeroth  moment  case  and  then 
relate  the  second  moment  evaluation  to  that  discussion. 

The  original  integral  has  the  form: 

»  f  ,  ,  ,  -C2t(z'  ‘  i£>T  1  <z'  -  jx)) 

‘o  ■  Cl/„’"  /.dz!  dz2'"  dzNe 


We  substitute  Z  =  (z  ‘  -  ju)  and  obtain: 


OO  CO 


*o  =  C1 


J  •••  /  dzidV‘*  dzNe 


-C2[ZX  R_1Z] 


.  CO  CO 


The  next  simplification  uses  the  fact  that  the  quadratic  form  is  unchanged  by 
a  coordinate  transformation  that  diagonalizes  the  matrix  R  \  This  transformation 

can  be  achieved  by  a  relation  involving  the  eigenvalues,  X  .  (i  =  1,  N)  and  corres- 

-1  1  T 

ponding  eigenvectors  of  R  .  Let  A  be  a  matrix  of  eigenvectors,  then  Y  =  A  Z 

•—  —  g 

reduces  the  quadratic  form  and  the  integral  to: 


Io  =  C,  /  •••  /  dyldy2 


dyNe 


-C2  tAlyl  +A2y2  +  ”  ’  ANyw'  ' 


For  convenience  we  transform  again  to  remove  the  eigenvalues  from  the  exponent: 


!0  =  C1  <XlX2  * '  •  V1/2  /•  •  ’/  d  W1  d  w2  * '  *  d  V 


-c2[wi2+w22+---wn2]1/2 


Up  to  this  point  the  relation  between  successive  sets  of  variables  has  been 
clear.  The  actual  evaluation  of  this  integral  requires  an  additional  coordinate 
transformation.  For  that  case,  the  Jacobian  will  have  to  be  considered  in  detail 
as  a  result  of  the  complexity  of  the  relations  between  the  two  sets  of  variables. 
This  is  discussed  in  Appendix  A  and  only  the  result  will  be  shown  here. 

8.  Hildebrand,  F.  B.  (1965)  Methods  of  Applied  Mathematics,  Prentice-Hall. 


'o  ■  C1(X.X2 


v-l/2/"rN-le-Srdrj'd9  / 

o 


•  N-2  . 

1  j  sm  Vld9N-l 
o  o 


/  sin  02  d02  or 


!0  =  (^1^2  V'1/2  27rC1C2N(N-l)<  /  sinN_Z  flN.j  dflN_j  /  01„  „  N_2 


sinN'3  0 , 


d0N-2'**  /  sin02d02* 


After  considerable  mathematical  manipulation,  a  closed  form  solution  is  obtain¬ 
able.  The  details  of  its  evaluation  are  presented  in  Appendix  B.  The  end  result  is 


I 

o 


=  2  (IT) 


N/2 


<V2 


r(N)Cj 


r(N/2>  =  i . 


(8) 


This  gives  us  the  first  of  the  two  relations  for  the  two  unknown  quantities  and  C2- 
For  the  other  relation  we  consider  the  second  moment  integral  for  each  of  the 
individual  variates: 

,  2  r  r  ,  2  -C2[(z'-H)T  R'1  (x'-H)]1/2 

I2(z')  =  ff  =  C1  ^  *  4  <v  Mj)  e . dz J  dz'  .  •  •  dz.jJ  . 

N 

Note  that  Ncr2  =  1  (z.  ).  After  a  change  of  variables  we  have: 


i=  1 


no2=  c.  f  •••  f  <z  z  +  z  z  +  .-*  +  z  2)e 

A-00  -  oo 


2_-c2!5t  51 


1/2 


dZ  •  • »  dZ  . 
1  N 


As  in  the  previous  case  we  next  decorrelate  the  variables  using  the  eigenvector 
matrix  transformation.  Here  though,  there  is  the  additional  multiplying  term  which 

oo  o  nr 

has  to  be  considered.  Note  that  +  z2  +  *  **  +  ZN  =  ZrZ  and  Z  =  A  jr  . 


Then 


=  (iT  4t>  =  yTy =  *  >r22  +  ***  +  yN2 

2  2 

If  we  again  use  the  relation  Wj  =  ^  we  have: 


No 


-c  as- 

e-c2[Wi 


«/W,2  w2 

+  —4-  + 


w. 


N 


-C9[W12  +  w22  +  •••  wN2] 


•  ■  •  4-  — — 

1/2 


dWj  ••• 


dW 


N 
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As  before,  the  relations  discussed  in  Appendix  A  are  used  to  carry  out  a  final 
transformation.  The  corresponding  detailed  evaluation  of  the  integral  is  covered 
in  Appendix  B  and  the  end  result  is 


2(,)N/2(X 


iV“  V'1/2r<N  +  2> S  V1  c1c2"(N'r2)/r(N/2). 

i=i  / 


This  is  the  required  second  relation  for  the  two  unknowns. 

N 

To  simplify  we  note  that  A.~*  is  the  trace  of  the  inverse  of  R~^  or 
N  i=l  1 

X.'J  =  Tr  (R-1)-1  =  Tr  (R)  =  NO2. 
i=l  1 

Substitution  and  rearrangement  leads  to  the  result: 


(2)  “2  (2ti) 


f  <n+i)n'2. 


C2  =  (N+ 1);  and  finally 

PE  <Z1*  z2*  ZN*  =  D 


1 1/2  N+  1 
R|  '  (2)  ~~T~  (2n) 


>  .. \T  d“  1  /„ *  .. 


v/2  -[(N+1)(*#-Mr  (10) 

(n+  i)  7  e  ~  ~  ~  ~  ~ 

We  have  now  arrived  at  the  forfn  of  our  N  -  variate  exponential  PDF.  Some 
final  comment  as  to  PE(z z'  . . . ,  z^)  satisfying  the  requirements  of  a  probability 
density  function  may  be  appropriate.  First,  the  zeroth  moment  relation  was  chosen 
to  satisfy  the  requirement  that  the  cumulative  distribution  over  all  space  be  unity. 
Secondly,  since  the  quadratic  form  is  non -negative  and  R  is  positive  definite,  this 
implies  0  £  PE(z',  z',  ....  z^.)  -  1.  Thus,  the  criteria  for  a  probability  density 
are  indeed  satisfied  by  this  form. 

The  determination  ot  the  appropriate  expressions  for  the  two  probability 
densities  in  which  we  are  interested  allows  us  to  proceed  to  the  actual  incorporation 
of  these  results  into  a  test  for  their  degree -of -fit  to  the  data.  This  will  be  the 
subject  of  the  next  section. 

4.  HYPOTHESIS  TEST 

As  has  been  discussed  in  the  introduction,  the  concern  for  our  particular  case 
is  to  determine  whether  the  various  terrain  subregions  are  better  described  by  a 
Gaussian  or  exponential  PDF.  This  binary  decision  could  be  generalized  depending 
on  the  data  involved  and  indeed,  could  even  be  made  a  more  complex  decision  pro¬ 
cess  if  the  appropriate  criteria  are  available. 


II 
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The  form  ox  the  hypothesis  test  used  here  is  based  on  the  maximum  a  posteriori 
4 

probability  criterion.  This  is  equivalent  to  a  minimum  error  probability  criterion. 
We  assign  hypothesis  H ^  to  the  Gaussian  case  and  hypothesis  lo  the  exponential. 
Then  the  likelihood  ratio  parameter. 


p_  (Z 
.A  1 

A-  nz 


-N' 


o 


1* 


Let  P  (H  )  be  the  probability  that  hypothesis  Ho  is  true.  Then  the  decision  rule 
may  be  written  as:  Choose  H.  if 

P(H  ) 

\  >  o 

A  ~  r^Fnrr  ■ 


For  our  case  we  assume  that  it  is  equally  likely  that  hypothesis  H,  or  H  is  true 

1  c 

and  the  decision  rule  reduces  to  whether  or  not  A  -  1.  Note  that  it  may  be  possible 
to  alter  the  probability  that  Hq  is  true  based  on  external  evidence  (such  as  the  type 
of  terrain). 

This  formulation  represents  the  decision  as  a  quotient  of  the  two  PDF's  of 
interest.  The  next  step  is  to  derive  the  specific  test  ratio  by  introducing  the  two 
expressions  for  the  N-variate  joint  PDF's  from  Eq.  (6)  and  Eq.  (10); 


N+l 

(2)  2  (2  «  ) 


N-l 

2 


N/ii  |  R| 1/2  (N+l) 


2  9  T  _  1 

where  again Q  =  (z  -  #x)A  fj  (z*-  Jf). 


After  some  manipulation  this  reduces  to: 


A 


P1  e(N+1)/2^  (Q  -  Vn+1)2 

r 


For  convenience,  we  rewrite  the  test  in  logarithmic  form  and  assert  that  Hj 
is  true  if  £n  X  -  0.  We  then  have  the  result:  Hj  is  true  if 


-  (Q  -  D2  -  \  £n  it  -  |  (N+  1)  +  ~  in 


r  mj.  i  1 


)  -  fn 


For  the  data  sets  used  in  the  example  of  this  report  N  =  100  so  we  can  translate 
the  decision  parameter  to  a  specific  value  for  this  case.  We  then  say  that  the 
terrain  heights  in  a  given  subregion  are  distributed  according  to  a  Gaussian  PDF 
provided 


9.  22  ^  Q  ^  10.  88 


or 


85.01  -  (z'  -  M)T  R'1  (z  -  m>  s  118.37 

and  correspondingly,  the  heights  are  better  described  by  the  exponential  PDF 
when Q“ >  1 18.  37  or  Q"  <  85.  01 

Since  the  hypothesis  test  has  been  formulated  in  terms  of  the  quadratic  form, 
Q“,  the  characterization  of  the  various  subregions  depends  on  the  determination  of 
that  quantity  for  the  respective  data  subsets.  The  next  section  addresses  this 
question  by  pointing  out  some  of  the  computational  complications  associated  with 
these  large,  jointly  distributed  data  sets. 


5.  POSITIVE  DEFINITENESS  OF  COVARIANCE  MATRIX  AND 

PROBLEMS  OF  INVERTING  LARGE  MATRICES 

It  is  important  in  these  studies  to  explicitly  demonstrate  the  positive  definite¬ 
ness  of  the  quadratic  form 

v  /  t  -1  v 

Q"  =  (z7  -  M)  R  (z  -  ju)  . 

9 

This  is  necessary  for  two  reasons:  (1)  The  positive  definiteness  of  Q”  insures  that 
the  functional  forms  for  both  the  Gaussian  and  exponential  N  -  variate  PDF's  are 
such  that  they  satisfy  the  requirements  of  a  true  mathematical  probability  density 
function,  (2)  the  quadratic  form  Q“  contains  the  inverse  of  the  covariance  matrix 
R  If  the  matrix  R  is  large  there  can  be  difficulties  inverting  it  on  a  digital 
computer.  If,  however,  R  is  both  positive  definite  and  symmetric  then  special, 
particularly  efficient  algoi  ithms  exist  for  the  inversion  process. 

We  will  first  address  this  latter  aspect.  For  the  purposes  of  the  proof  we 
consider  a  more  general  form  which  involves  the  matrix  R  rather  than  R-1; 


S2  = 


T  „  * 
u  R  u 
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j 

! 


where  u  is  any  N-dimensional  vector  and  (#)  denotes  complex  conjugate.  Further¬ 
more,  by  definition,  S2  is  positive  definite  if  u"^  R  u  -  0,  where  the  equality 
holds  if  and  only  if  u  is  identically  equal  to  the  zero  vector. 

Now,  as  the  first  step  in  showing  that  S  is  positive  definite,  we  use  the  well 
known  formula  for  the  Fourier  transform  of  a  Gaussian  function 


J  e  at  e^1  dt  =^i(n  la)  t. 


u>2/4cr 


(ID 


where  Re  (or)  >  0. 


2  -Tmn/T2 

This  is  introduced  as  a  result  of  our  assumption  that  Rmn  =  a  6 


Relating  the  Gaussian  form  to  our  correlation  case  we  note  that  the  distance 
between  the  ordinates  is  equally  spaced  in  the  x  -  y  plane  so 


(ym  -  yn>2  -  'Ay>2  (rn  -  n)2 


and  likewise  for  the  distance  between  abscissae 


(x  -  x  )2  =  (A  x)2  (m  -  n)2  . 
m  n 


Let  a1  =  (1/4)  (T/Ax)2  and  a2  -  (1/4)  (T/Ay)2  so  that 


-(x  -  x  )2/T2 

m  n 


_ _  .«  -  a.  t.  i(m  -n)t1 

■^(o-j/tt)  /  e  e  dtj 

«*  OO 


(12a) 


and 


-(y„  -  yJ2/T2 


m  n  =  V  («,/ 'sr)  /  e 


«o  ~a2t2Z  i(m  •  n)tj 


^tn 


(12b) 


o  T  * 

=  u  R  u 


N 

N 

II 

M 

E 

m=l 

n=l 

N 

N 

-E 

E 

m=l 

n=l 

(yfo 
=  1 

;l“2 

7t 

•  -<xm  -  V2/T2  -  V2^2 

u  u  e  e 

m  n 


i(m  -  nH1  i(m  -  nH2  "  ^2^2 

f  dt«  e  e  e 

-  CO 


m=l 


n=l 


(13) 
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This  Leads  to  the  result 


’  u 


“2  ° 


2 —  i  dti  J  dt2 


"l4!2  ~a2X2  I 


N 

E 

£=1 


i£  (t.+t9) 


Uj  e 


(14) 


o  2  .....  . 

and  thus  S  -  0.  It  remains  to  be  sltown  that  S  =  0  iff  u  =  0.  Trivially,  if 

N 

u  =  0,  S2  =  0.  Also,  if  S2  =  uT  R  u  =  0  then  ^  u^  e1^1  =  0  follows  from  above 

where  t  =  t,  +  t9.  Next,  introducing  successive  complex  conjugates  of  the  exponen- 
1  1  2n  N 

tial  and  integrating  yields  f  e’^1  e  ***  dt  =  2tt  0^.  Thus  ^2  =  0  implies 

o  2  £~1 

u  =  0  for  k  =  1,  2,  •  •  • .  N  so  that  S  =  0  iff  u  is  equal  to  the  zero  vector  and  we 

K  "  o 

have  completed  the  proof  that  the  quadratic  form  S  is  positive  definite.  (Part  of 

9 

this  proof  is  based  upon  the  work  of  Fante.  ) 

ft 

S  being  positive  definite  implies  that  the  eigenvalues  A.  of  the  matrix  R  are 

all  positive.  Also,  if  Aj  is  an  eigenvalue  of  11,  i=  1,  ....  N,  then  the  eigenvalues 

of  R-1  are  Af1  and  these  are  all  positive.  Thus  a  quadratic  form  associated  with 

R“1  will  be  positive  definite.8  Since  Q2  =  (z'  -  m)T  R_1  (z'  -  m)  is  an  example 
~  2  ~  ~  ~  ~ 
where  the  vectors  are  real,  then  Q  is  positive  definite. 

The  next  topic  is  that  of  matrix  inversion.  The  data  sets  analyzed  in  this  re¬ 
port  were  such  that  N  =  100,  so  that  the  matrix  R  to  be  inverted  was  of  order 
100  by  100  elements  having  an  extended  range  of  magnitude.  This  combination  can 
introduce  complications  to  the  inversion  process. 

The  first  technique  used  for  inversion  was  a  Cholesky  reduction  or  square  root 
method  (Forsythe  and  ?*Ioler,  pp  14  et  seq).  The  specific  algorithm  that  was 
used  required  that  R  be  symmetric  and  positive  definite.  Both  of  these  criteria 
are  satisfied.  Because  the  matrix  R  was  particularly  large,  round-off  error 
accumulated  during  the  Cholesky  reduction  was  so  large  that  the  verification  of  the 
inverse  based  on  R  •  R  1  =  1.  the  unit  matrix,  resulted  in  a  matrix  with  off  diagonal 
elements  sometimes  greater  than  one,  instead  of  10 9  10 *  or  less  which  would  be 
expected  for  the  computer  used. 

To  circumvent  this  loss  in  accuracy,  the  more  conventional  Gauss  elimination 
method  with  pivoting  was  used  (Forsythe  and  Moler18).  Also,  double-precision 
accuracy  was  retained  in  the  computational  procedure.  When  this  also  had  unreli- 

_3 

able  results  in  some  cases,  a  threshold  value  of  10  was  set  so  that  any  element 
in  the  R  matrix  less  than  or  equal  to  the  threshold  was  automatically  set  equal  to 
zero.  When  these  procedures  were  applied  to  the  inversion  problem,  R  times  R  1 

9.  Fante,  R.  L. ,  Private  communication. 

10.  Forsythe,  G.,  and  Mole r,  C.  B.  (1967)  Computer  Solution  of  Linear  Algebraic 

Systems,  Prentice-Hall,  pp  114  et  seq. 
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yielded  a  matrix  whose  diagonal  elements  were  equal  to  1  ±  e,  where  t  ~  10 
and  the  off-diagonal  elements  were  ~e. 

There  are  numerous  techniques  for  improving  the  accuracy  when  obtaining  the 
inverse  of  a  matrix.  References  Forsythe  and  Moler1®  and  Berezin  and  Zhidkov11 
discuss  at  least  a  dozen  methods.  The  choice  of  a  particular  method  depends  upon 
the  peculiar  properties  of  the  original  matrix  (that  is,  the  order  of  the  matrix, 
positive  definiteness,  symmetry,  size  of  matrix  elements,  and  so  on).  Another 
factor  in  the  choice  of  an  inversion  method  depends  upon  trade-offs  in  computer 
storage  requirements  vs  computation  time. 

In  applying  the  procedure  outlined  in  this  report  to  other  data  sets  and  PDF's, 
the  appropriate  forms  of  the  correlation  of  the  variates  have  to  be  considered  and 
individual  decisions  made  as  to  inversion  procedures  where  necessary. 


6.  STATISTICAL  CHARACTERIZATION  OF  A  SPECIFIC  SITE 

The  specific  example  used  in  this  report  is  a  site  in  Eastern  Massachusetts 

12 

selected  because  electromagnetic  scattering  data  from  the  terrain  is  available. 

A  rectangular  area  around  the  Discrete  Address  Beacon  System  (DABS)  site  was 
designated.  The  rectangular  area  w«s  43.  3  km  long  and  42.  8  km  wide.  The 
rectangular  area  was  then  subdivided  into  smaller  rectangular  cells,  each  with 
sides  of  2050  m  by  1852  m.  Each  cell  is  then  further  subdivided  into  a  10  by  10 
grid  of  points.  A  data  base  of  topographic  elevations  for  this  area  is  available  at 
the  Electromagnetic  Compatibility  Analysis  Center  (ECAC).  This  was  prepared 
from  Defense  Mapping  Agency  (DMA)  supplied  digitized  terrain  maps  at  1:250,000 
scale  size. 

The  statistical  analysis  techniques  discussed  in  Sections  2  and  4  of  this  report 
were  then  applied  to  each  cell,  where  now  (z/j  represents  the  set  of  topographic 
elevations  at  the  grid  points  in  each  cell  (i  =  1,  2,  3,  ....  100).  In  addition  to 
parameter  estimates,  each  cell  has  a  geologic  feature  code  associated  with  it. 

Table  1  indicates  the  particular  geologic  codes  used. 

Table  2  indicates  the  results  of  applying  the  statistical  procedures  indicated  in 
this  report  to  the  DABS  radar  site.  The  rectangular  region  was  subdivided  into 
528  cells,  with  twenty-two  cells  extending  along  the  x-axis  of  the  grid  and  twenty-four 
cells  extending  along  the  y-axis.  The  (x,  y)  coordinates  of  a  particular  cell  are 
referenced  with  respect  to  the  center  of  that  cell.  The  origin  of  the  (x,  y)  coor¬ 
dinate  system  is  taken  at  the  center  of  the  cell  in  the  extreme  southwest  corner 

11.  Berezin,  I,  S. ,  and  Zhidkov,  N.  P.  (1965)  Computing  Methods,  Vol.  II, 

Pergamon  Press. 

12,  McGarty,  T.  P.  (1974)  Models  of  Multipath  Propagation  Effects  in  a  Ground-to- 

Air  Surveillance  System,  Ab  777241. 
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of  the  rectangular  region  (see  Figure  1).  In  Table  2,  the  columns  labelled  X  and  Y 
refer  to  the  (x,  y)  coordinates  of  a  particular  cell.  The  indexing  of  the  cells  is 
such  that  as  one  moves  down  the  X  and  Y  columns,  increments  equal  to  one  cell 
length  in  the  x-direction  are  taken  with  the  y-eoordinate  fixed  until  the  extreme 
eastern  edge  of  the  rectangular  region  is  reached  (see  Figure  1).  Then,  the 
sequencing  returns  to  the  extreme  western  edge  with  the  y-coordinate  increased  by 
one  cell  width  and  increments  again  are  taken  in  the  x-direetion.  This  procedure 
is  repeated  until  the  extreme  northeast  corner  is  reached. 

The  column  labelled  1C  contains  the  geologic  codes.  The  column  labelled 
SIG  lists  the  mean  heights.  VAR  refers  to  the  variance  in  surface  height,  TA  to 
the  correlation  length  and  TEST  to  the  results  of  the  hypothesis  testing  procedure. 
The  units  of  length  are  in  meters. 

It  may  be  recalled  from  Section  4  that  (since  TEST  =  Q2  -  118.  3?)  the  heights 
in  a  given  region  are  Gaussian  distributed  when  -33.  36  £  TEST  -  0. 0  and  exponen¬ 
tially  distributed  otherwise.  When  TEST  =  0.  0  and  VAR  is  very  small,  the  region 
is  essentially  a  smooth  surface  (no  roughness). 

One  observation  that  can  be  made  is  that  when  the  magnitude  of  TEST  is  very 
large,  so  is  the  correlation  length;  TA  for  those  cases  is  comparable  to  one-half 
the  cell  size  or  even  larger.  When  that  occurs,  the  determinant  of  the  covariance 
matrix  R  becomes  very  small.  As  a  result  it  becomes  more  difficult  to  obtain  an 
accurate  inverse  of  R  due  to  rapid  build-up  of  round-off  error.  A  related  problem 
with  the  results  is  that  for  some  cases  TEST  <  -118.37.  This  could  only  occur  if 
the  quadratic  form,  Q2  were  not  positive  definite.  This  contradiction  of  that 
theoretically  imposed  condition  suggests  the  presence  of  additional  machine-induced 
errors.  Those  results  can  not  be  valid.  Further  investigations  of  these  aspects 
are  being  conducted,  along  with  assessment  or  the  parameter  estimation  techniques. 

In  order  to  use  these  results  in  the  rough  surface  electromagnetic  calculations, 
one  additional  aspect  should  be  noted.  For  the  types  of  geological  features  that 
describe  the  respective  regions,  data  exists  on  the  associated  complex  dielectric 
constants  at  microwave  frequencies. ^ 


13.  Lytle,  R.  J.  (1974)  Measurement  of  earth  medium  electrical  characteristics: 

Techniques,  results  and  applications,  IEEE  Trans,  on  Geoscience 
Electronics,  GE.- 1,2:81. 

14.  Long,  M,  W.  (1975)  Radar  Reflectivity  of  Land  and  Sea,  Lexington  Books. 
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Table  1.  Geologic  Code  for  the  Site 


0  -  Sea  Water 

4  -  Cleared  Area 

1  -  Woods 

5  -  Town  (City) 

2  -  Swamp 

G  -  Road 

3  -  Lake 

7  -  Intermittent  Lake" 

The  Massachusetts  area  has  few  true  inter¬ 
mittent  lakes.  The  code  was  used  when  many 
small  lakes  or  ponds  appear  in  a  cell  (rather 
than  a  large,  all  water  cell). 


RtCTASGULAR  REGION 


Figure  1.  Subdivision  of  Site 
Into  Sample  Regions 
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Table  2.  Details  of  the  Statistical  Analysis  (Cont) 
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Table  2,  Details  of  the  Statistical  Analysis  (Cont) 


Y  _  IC  SIG  VAR _  TA  TEST 


o.a 

11109.6 

7 

,  248  '4-03 

. 8556+03 

,  5  2*8  5+0  3 

-.7596+05 

2051  . i 

11109.6 

•4 

.l*?***? 

. H34E+03 

.648E+03 

.115E+06 

4102. 2 

11109.6 

K 

•  152E+P3 

.*50E+03 

.9191+03 

6157.3 

11109.6 

y 

•  1  06r+0? 

.2905+03 

.6405+03 

•  220E+06 

5704,4 

11109,6 

1 

. 1 06*  +  07 

.18  75+03 

.5ieE+03 

‘.'3196+04 

10255.5 

11109.6 

1 

.1175+0? 

,9285+02 

.531E+03 

-.2726+04 

12306.6 

11109,6 

u 

.8015+0? 

•  1885  +  03 

.5605+03 

-  .2301+05 

14353.? 

11109.6 

5 

.871  +  4.02 

,5495+03 

. 688E+03 

.1606+06 

icico.e 

11109.6 

.111**-'? 

.1155+03 

.3255+03 

.2446+03 

13459.5 

11109.6 

u 

« 127-+'  * 

.  9f  16  +  1.? 

. 558E+03 

.68  -6+04 

20511.0 

11109.6 

■4 

.1045+07 

. 8716+03 

.739E+03 

.2836+07 

22562.1 

11109.6 

u 

.  o  ■»«,<■  +  07 

.1986+03 

•  460E+Q  3 

.6375+03 

24613.2 

11109.6 

7 

.76»*+0? 

.3125+02 

.435E+03 

.106E+Q4 

25664.3 

11109.6 

c 

.?6*r+J2 

.  5755  +  02 

.352E+03 

.3495+03 

20715.4 

11109.6 

i 

. 74*5+02 

.5475+02 

•  4  67E+0  3 

.7055+04 

30766.5 

11109.6 

i 

.  146E  +  03 

.427E+03 

-.438E+04 

32*17.6 

1 1109. 6 

i 

.•6465+02 

. 922E+E2 

.4805+03 

.1545+04 

3^863.7 

11109.6 

i 

4-07 

.1166+0? 

.6735+0? 

.1795+06 

35915.8 

11109,6 

i 

. 611 *  +  02 

, 55  4£+  01 

.  +6  2E+Q  3 

-.7865+03 

33970.5 

11109.6 

4 

.■8?r+52 

. 4?  7E  +  Q2 

.6875+03 

-.635E+06 

41022.0 

11139.6 

5 

*  ^  27  r+ r  ? 

•  <1 1E+C2 

.5456+03 

.3375+84 

43073.1 

11109.6 

5 

.4646+0? 

.  39  36  +  02 

•520E+03 

.1505+04 

6.0 

12961.2 

7 

. ?316f()7 

.  2255+04 

.6765+03 

-.685E+07- 

2051. i 

12961.2 

1 

.1356+07 

.178E+03 

.4855+03 

.5295+04 

4102.2 

12961.7 

C 

.  1  005+  O'* 

.6655+02 

•372E+03 

.4436+03 

5153.3 

12961.2 

* 

.127-+03 

.  20  85  +  01 

.441E+03 

.786E+0  3 

3204.4 

12961.2 

7 

♦  ii** 

.8125+02 

.7765+03 

.2075+07 

10255.5 

12961.2 

7 

.1  14; +'-3 

. 5086+02 

.4485+03 

.552E+04 

12306.6 

12961.7 

\ 

.7666+02 

.4025*02 

.5845 +03 

.1545+05 

14357.7 

12961.2 

,9  11 '*  +  0? 

.157E+03 

.6425+03 

.2265+06 

15400.8 

1 2°61 , 2 

*7 

.7935+0? 

.2266+03 

.4545+03 

.4^55+0  4 

13455.9 

12961.2 

5 

,8966+0'> 

•  11 *E+03 

.468E+03 

.5935+04 

23511.0 

12961.2 

1 

.1065+03 

.6115+03 

•420E+03 

.4725+04 

22552.1 

12961.2 

1 

•7Ti Hf r? 

. I7  TE+C3 

.359E+03 

.3515+03 

24613.2 

12961.2 

l 

.  056J+6? 

.  E-4E+C2 

.368E+Q3 

.412E+03 

26664.3 

12961.2 

u 

.841 6+82 

.1366+03 

.3345+03 

•603E+03 

23715.4 

12961.2 

1 

,69r*+07 

.2196+02 

• 467E+03 

.7415+84' 

30766.5 

12961.2 

1 

. 6?7r+  02 

.112E+03 

. 635E+03 

•361E+07 

32817.6 

12961.2 

7 

.6775+62 

, 66 1E+C2 

.5215+03 

.313E+05" 

34363.7 

12961.2 

7 

.6566+02 

.68  56  +  02 

.5345+03 

.663E+04 

36515.8 

12961.2 

l 

.514E+6? 

.  467E+02 

•567E+03 

-.205E+05 

33970.9 

12961.2 

7 

. 4035+02 

,  222E+02 

. 4345+03 

-.508E+03 

41072. 6 

12961.2 

- 

,4376+0? 

. 164E+02 

• ?75E+D3 

.634E+03 

43073.1 

12961,2 

s 

.*445+0? 

.  50  76+02 

.569E+03 

-.562E+04 

3.G 

14812.8 

c 

.2066+0? 

.1036+04 

•725E+03 

.2895+07 

2051.1 

14812.  8 

- 

. i’g?  +  n3 

■  3166+03 

.4676+03 

•347E+04 

4102.2 

14812.8 

c 

.1148+0? 

. 2766+03 

.5525+03 

.3515+04 

5153.3 

14812,8 

5 

.1406+03 

.  39 5E+03 

. 664E+03 

.30  95+06 

3204.4 

14812.8 

7 

.1056+07 

.3975+02 

.5705+03 

.2485+05 

.  10255,5 

14812.8 

7 

#cccc*.n2 

. 5746+02 

.4675+03 

.8335+04 

12305.6 

14812.8 

l 

.7525+02 

,9236+02 

. 429E+03 

.7915+03 

14357.7 

14612. 8 

l 

.1186+03 

.554E+03 

.1355+03 

16403.8 

14812,8 

i 

,7646+82 

. 2936+02 

.4526+03 

.8355+04 

18459.9 

14812.3 

.8606+02 

.34  45+02 

.6305+03 

.6575+06 

2051  i.O 

14812,  8 

7 

.9006+0? 

,0176+02 

.435E+03 

.22HI+0Y 

22562.1 

14812.8 

7 

.78*6+02 

.  20  85  +  03 

.4525+03 

.8885+03 

2461 3".  2 

14812.8 

A 

.72«c+r? 

.5535+02 

.4445+03 

•768E+04 

26664.  3 

14812.8 

7 

.*395+0? 

.8265+02 

. 356E+03 

.5175+03 

23715.4 

14812,8 

I 

.6446+1;? 

.  5975+02 

.565E+03 

-.3385+04 

33766.5 

14812.8 

7 

.6255+02 

, 45 15+ C 2 

. 540E+03 

.1375+05 

32817.6 

14812. 8 

7 

.6745+0? 

.6055+02 

.51  gE+OlS 

-  .  48TT+1T 

34863,7 

14812,3 

4 

.607--,  02 

.61 6E+02 

.544F+03 

.292E+05 

36919.8 

14812.8 

5 

.5715  +  02 

.  38  9E+02 

.5275+03 

-.8395+05 

38970.9 

14812,  8 

5 

,5425+02 

.  624E+02 

.7255+03 

.3075+07 

41022.0 

14812.8 

5 

.4776+0? 

.118E+03 

.4225+03 

-.1525+8?' 

43073,1 

14812.6 

1 

.6045+02 

.479E+C2 

. 434E+03 

-.103E+04 
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Table  2.  Details  of  the  Statistical  Analysis  (Cont) 


X 

Y 

IC 

SIG 

VAR 

TA 

TEST 

0.0 

1  6664. 4 

5 

.17154-03 

.4275+07 

.873E+03 

-.7626+08 

’051.1 

16664.4 

c 

.1176303 

.  ?9?E+03 

.461E+03 

.7496+04 

4102.2 

1  6664*4 

5 

,  1*7  r  +  03 

.2696*03 

, 376E+03 

, 934E+0  3 

6153. 3 

1 6664_.  4 

c 

.1716303 

.1616+03 

.521E+03 

.3416+03 

0204.4 

1.6664',  4 

1 

. 96°"3  *2 

.  27  5E+C2 

.  447E  *03 

-.4286+04 

13255 *  5 

16664,4 

1 

.1016303 

.  3976+02 

.5416+03 

•234E+05 

12306.  6 

16664,4 

1 

.  23  °E+  C-3 

.8066+03 

.II6EV08 

14357,7 

16664.4 

1 

.724-30? 

.  83  8F  +  01 

. 341E+Q  3 

.4136+03 

16408.8 

16664 ,4 

1 

.71*8  63  02 

.  20  ?E  +  02 

. 741E+03 

.3486+03 

10459.9 

1_6664*  4 

1 

*77t;=>g? 

.677E+02 

, *976+07 

.131E+04 

23511  * C 

16664. u 

7 

.  77" F*  ?  2 

.  51C5  +  C-2 

.4186+03 

.2506+04 

22562  •  1 

1 6664*4 

1 

. 7?7r+c2 

•  1+  86+02 

.4326+03 

-.8086+03 

24613.2 

16664.4 

y 

,7-J7Cfl)2 

.1756+02 

.361E+Q3 

-.9546+02 

?666<i  *3 

16664,4 

? 

.728730? 

.28 5E+02 

.5236+03 

.173E+05 

2fl/l 5*4 

16664.4 

1 

.7°?F+02 

.  22  9-6  +  0  3 

. 52  5E+03 

.342E+05 

33766*  5 

^6664*  4 

1 

.■’1'>73"c 

.6176  +  111 

.560E+03 

-  .273E+05 

32817,6 

16664.4 

? 

.6496302 

.  4136  +  02 

.3846+03 

.757E+04 

348.60.7. 

16664.4 

c, 

*604^4 *2 

.2276+02 

.3826+03 

.347E+Q4 

36919.8 

16664.4 

5 

.5486*02 

•  398E+C-2 

. 4  33E+03 

.872E+03 

39170.9 

16664. t, 

c 

,f ?7f *02 

.  44  7E  +  02 

.3656+03 

.5886+03 

41022,0 

16664,4 

5 

, 440730? 

.8206+02 

.  787E+03 

•  343E+0  7 

wozia 

1£6M,4 

1 

.  61°  r3f2 

.8456+02 

•53SE+03 

.1552+1& 

3*0 

18516.  C 

7 

.15473"? 

. 5436+03 

.  6476+03 

.1076+06 

2051.1 

18516.0 

.1256303 

.1656+03 

. 654E+03 

-.1656+07 

4102.2 

18516.0 

1 

•  15’ 630 7 

.707E+02 

. 466E+03 

•172E+04 

.  £15.3, .? 

18516.0 

1 

.118630? 

. 274E+03 

.53?E+03 

-.6336+02. 

8204,4 

18516.0 

1 

.110630-* 

.1666+03 

.5446+03 

•661E+04 

iM5.5t5 

185.16.0 

1 

,107730* 

.8266+02 

.3436+03 

.328E+03 

1233  *6 

18516. C 

1 

, cp7Cf r? 

. 143E+C3 

.5916+03 

.4046+05 

14357,7. 

19516.0 

1 

.7i?7,<-£ 

.  631E+C2 

. 560E+03 

-.3196+04 

1640  *!* 

18516.0 

1 

.?4F  7*0? 

.4526+02 

•760E+03 

.474E+03 

i.9  459 , 9 

185.16,0 

4 

.7557*02 

.295E+03 

•438E+03 

•550E+03 

20511.0 

18516.0 

1 

.9097*02 

.1816*03 

.4246+03 

-•381E+05 

22512,1. 

18516,0 

7 

.767  7*0? 

•F06E+02 

.5686+03 

-  .1346+06 

2461 3*2 

18516.0 

1 

.6877*02 

.3826+02 

,4326+03 

-.11ZE+04 

26554.3 

18516.0 

1 

•677c*r? 

, 17PE+02 

.462E+03 

•  114E+0  4 

28715.4 

16516.0 

1 

•  977  r*  62 

.9846+02 

. 457E+03 

•98EE+03 

.10766.5 

18516.0 

5 

.7367*419 

.1615+02 

.416E+03 

.199E+04 

3781  7  ,6 

18516,0 

1 

,5836*0? 

.1516+02 

.  434E+03 

•  224E+0  3 

34868, 7 

18516,0 

c 

.6246*02 

.60 1E+02 

. 322E+03 

-•354E+03 

36919,8 

18516. n 

5 

.  44263  r  2 

•  1966+C  2 

.49  0E+03 

,4776+04 

38970,9 

18516.0 

5 

,5657*02 

.  6236+02 

•437E+03 

-.1866+03 

41022,0 

16516. C 

1 

, 776  7*  n? 

.  533E+C-H 

.643E+03 

•56ZE+05 

43073.1 

18516,0 

1 

.4  447*0? 

.  P03E+02 

.  635E+03 

.537E+06 

0.0 

70367,6 

7 

.1586*03 

. 700E+03 

•543E+03 

•126E+05 

.  205.1a 

?.036?i6 

1 

,1*87*03 

,  23  4E+03 

,3706+03 

.507E+Q2 

4102.2 

20367.6 

1 

,1667+r* 

.104E+03 

•443E+03 

.723E+0 3 

-615.3.3 

20367.6 

1 

.54°"+"? 

.  2436+03 

.6756+03 

•  ?Z5E+Jit 

8204.4 

20367.6 

7 

•I?’"*"* 

,1516+03 

.4805+03 

•1O1E+04 

10255,5 

20367.6 

1 

,  31 56  +  03 

.7706+0? 

-»6l4t+07 

12306.6 

20367.6 

1 

.107630? 

,  661E+02 

. 441E+03 

•24EE+04 

14757.7 

20367.6 

1 

.’13730? 

, 368E+0? 

.350E+O7 

•426E+03 

16408.8 

20367,6 

1 

.7857+02 

.  206E  +  03 

•  7  066+0  3 

.1116+07 

18459.9 

20367.6 

5 

.1016307 

.  201E+03 

•453E+03 

20511.0 

20367.6 

1 

, 9777+0? 

•  28  ?E  +  C  3 

. 657E+03 

-.6016+06 

22562.1 

70367,6 

1 

• '^4c ♦ 0? 

.592E+02 

•481E+03 

-.6016+04 

24613,2 

20367,6 

i 

.7827+0? 

.5216+02 

•610E+03 

.6456+05 

26664.3 

£5.361,5 

1. 

.6766*02 

.  2066+03 

• 42iE+2 3 

*1  +  76+0  4 

28715,4 

20367.6 

c 

.7395*02 

.2966+03 

•647E+03 

.2746+06 

10,76.6, ,5 

2.0367.6 

5 

.660F30? 

,  8506  +  02 

.’53E+03 

.344E+CL3 

32817.6 

20367.6' 

c 

.7246*02 

.  164E  +  C3 

.457E+03 

.280E+03 

34868.7 

2  0367.  6 

5 

..6P6  7+0? 

. 157E+03 

• 397E+03 

.7732+03 

36919.8 

20367.6 

I 

,1907*0? 

.  884E  +  01 

.6566+03 

.4636+05 

38970.9 

20367,6 

1 

.474  -  +  n? 

,  1286  +  02 

.615E+03 

•  161E+C6 

41022.0 

20367,6 

1 

.3226*02 

.134E+02 

. 569E+03 

-.7466+04 

43C73.1 

20367.6 

7 

.3967+7? 

,4946+02 

•599E+03 

•102E+06 
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Table  2,  Details  of  the  Statistical  Analysis  (Cont) 


X 


0.0 

2061.1 
4102.2 
,615?.’ 
8204  •  4 

10255 . 5 

12306.6 

14357.7 

16408.8 

18459.9 
20511.0 

11562.1 

24613.2 

26664.3 

28715.4 

30766.6 

32817.6 

34868.7 

36919.8 

33973.9 

41022.0 

43073.1 

0.0 

2051.1 

4102.2 

6153.3 

3204.4 

10255.5 


17306.6" 

14357.7 

16408.8 
19469.  g 
205 11,0 
11562.1 
2461 Y.2 

26664.3 

23715.4 

33766.5 

32817.6 

34868.7 


36919.8 

39970.9 
41022.3 

43073.1 
3.0 

,2051.1 

4102.7 

6153.7 
3204.4 

13255.5 

12306.6 
A.4757.7 
1 64  0*9  •  8 

19459.9 
20511,0 

22562.1 

24613.2 
,26664.7 
23715“.  4 
30766,5 
32°  17. 6 
34863  .7 

36919.8 

39973.9 
41022.0 
43073,1 


y 

IC 

SIG 

VAR 

TA 

22219.2 

T 

.186=4-0? 

. 126E+04 

.5286+03 

22219.2 

1 

.16=  E+03 

.2?7E+C3 

.6526+03 

22219.2 

1 

.  1=0  =  4'.? 

.91FE+12 

.4516+03 

72219.2 

1 

.l?i=+P? 

•368E403 

.  652E+0  3 

22219.2 

7 

.117=407 

, Q?lE+02 

.534E+03 

22219.2 

1 

.ll°F+03 

.2256+03 

.5506+03 

2 22lS. 2 

l 

.Q2==+r.  2 

.27=6+63 

•633E+03 

22219.2 

l 

.7745*0? 

.74  0E+02 

.439E+03 

72219.2 

i 

.?0“F+02 

. 111F+C3 

•  7746+03 

22219.2 

.1  07=403 

. 101E+03 

. 675E+03 

72210.7 

1 

.antE+n? 

.  355E  +  03 

.761E+03 

22219.2 

1 

. 707=40? 

, 12?£+03 

.547E+03 

22219.2 

1 

.814=4=2 

.  2?CE*  03 

•  346E+0  3 ‘ 

22219,2 

7 

.S8=c*:2 

.  25  CE  +  C3 

.438E+03 

22219.7 

7 

.595=40? 

•  67  6E  +  01 

.672E+03 

22219,2 

5 

,576=402 

.  132E+03 

•367E+03 

22219.2 

5 

.707=402 

.  38  7E+03 

.5186+03 

22219.2 

c 

.6615+02 

•196E+03 

.576E+03 

22219.2 

f 

.303=40? 

.  916E  +  01 

• 636E+03 

22219.2 

c 

.741=4=2 

. 15  6E+C2 

.5256+03 

22219.  2 

i 

.’4:=+C5 

.  292E  +  C2 

.5576+03 

22219 . 2 

t 

.467=40? 

.  757E  +  02 

•4456+03 

24070,6 

1 

.?4<r+0? 

•894E+03 

. 644E+0  3 

24070.8 

1 

.187=407 

.54 0E+03 

•629E+03 

24070.8 

1 

« 1 =7=4  =3 

. 296E+03 

.7396+03 

24070.8 

1 

.132=403 

.1P1E+03 

. 594E+03 

74070.  8 

1 

.il°  =  4C7 

•153E+C3 

•451E+03 

24070,8 

1 

.107=403 

.  28fE+C  3 

.572E+03 

74070.8 

« 

,^ec c+n? 

.?9  PE+03 

•615E+D3 

24070.8 

t 

|H(s=4.0? 

.3266+03 

, 770E+03 

24070.8 

7 

.53-  =*02 

. 149E+02 

.5836+03 

24070.8 

7 

.047=4*2 

.15 1E+03 

.6476+03 

24070.8 

1 

•6?6=+Q2 

.2=5E+0? 

.3626+03 

74070.8 

1 

.7=4=4=? 

.  81 5E+C2 

.4756+03 

24070.8 

7 

,»78f+02 

.  39  6E+D3 

.  6646+03 

74070.8 

1 

.771=40? 

« 10  2E  +  03 

.471E+03 

24070,8 

1 

.f67=*02 

,  370E  +  02 

.4296+03 

74070.8 

5 

.5==  •*•’2 

•  253E+02 

.3756+03 

2407C.  8 

c 

,419=4'? 

.878E+61 

.3416+03 

74Q70.8 

c 

• 41cc*o2 

.  592E+02 

.4506+03 

74070.8 

5 

.415=40? 

.104E+02 

. 524E+0  3 

24070.  8 

*5 

.329=402 

.  393E+02 

. 642E+03 

24070.8 

.772E40'» 

.  8356+0? 

•432E+03 

24070.8 

7 

.42?=*0? 

.18  06+0? 

.5296+03 

25922,4 

7 

.224=48? 

.636E+0? 

.5876+03 

2  2  c  %  u 

1 

.174=4=3 

.  65CE+6 3 

.4546+03 

25927.4 

1 

.179=407 

.  28  36  +  03 

. 779E+03 

25922.4 

1 

.1?==407 

.  316E  +  03 

.5226+03 

25922,4 

1 

. 1075*07 

.  71CE+03 

.6356+03 

25922.4 

7 

.=616+02 

.5256+03 

25922.4 

1 

,104=*03 

•  160E+03 

•561E+03 

25922.4 

i 

.=34=4=7 

.  278E+C3 

.6206+03 

25922.4 

1 

.68?i*r2 

.157E+8? 

. 576E+03 

25927.4 

7 

.7=7r  4112 

.  115E  +  03 

.6356+03 

2  59  22,4 

7 

. '34=40? 

.  25=6+03 

.4736+03 

25922.4 

1 

.8  05=402 

.8126+02 

.  7546+03 

25922.4 

7 

.  77-  r  .r  3 

•255F+03 

.5246+03 

25922.4 

l 

. 7ot=. 02 

.116E+03 

.  7526+03 

26122.4 

9 

,=6*54"? 

.5166+02 

.4866+03 

TEST 


26922.4 

75922.4 

25922.4 

25922.4 

25922.4 
2=922.4 

25922.4 


•  586 '  +  02 

, 737r*P7 

,?94=+0? 

.•»07?*02 

•268T+D2 

,700=402 

•  4  =8  =4 ' 5 


. 176E+13 
.627E+0? 
.  4°  76+02 
•?,cE+02 
.  52  86+62 
.  1476+03 
.123E+C2 


.2626+05 

.1D7E+06 


•1066+04 
•  40SE+06 
.2176+04 
*•2466+05 
-.6236+06 
.??7E+0_3 
•307E+07 
.333E+06 


. 6286+03 
.5486+03 
.  613E+03 
.5386+03 
.718E+03 
.45  IE  <-0  3 
•  3786+03 


.742E+07 

-.627E+05 

•3626+03 

•7846+03 


•1036+0  6 
.5336+03 
.727E+0  4 
»15SE+0J^ 
>. 2911+06 
•138E+05 
-•209E+04 
•9556+04 
•521E+05 
.2756+06 
“.1166+07 
.229E+06 
-“.1826+04' 
.4996+05 
,  .1596+06 
•458E+07 

.z?ir+04 

•7996+05 

•529E+0T 

.617E+03 

•416E+06 

.136E+05 

-»72lfe+03 

.721E+03 

.4116+03" 

.196E+04. 

.121E+D5 

♦171E+06 

.115E+04 

-.900E+04 

.675E+0~5 

•141E+04 


.357E+07 

.128E+05 

.1616+07 

.2046+05 

.5746+04 

•  339E+06. 
.1536+05 
.1696+07 

-.229E+04 

.195E+07 

.9696+04" 

,5936+06 

•  848E+"0  3_ 
.6006+0,5 

-.114E+05 

•  655E+0  5. 
•252E+05 

•  29  2E+07 
•  •  1 0  5E  +0*4 

.4286+03 


?TKSFaG2IS 

JiwM  'COz*  x  -  '■'* 
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Table  2.  Details  of  the  Statistical  Analysis  (Cont) 


o.o 

’=?  1.1 

61<H.’ 
’2=4.  b 
10’55.5 
l’3fi6.fc 

147=;’. 7 
164' 0, 8 
1 ’455.0 
23=11.0 
2’==2,1 
’461 ’.  2 
’6664,  7 
207i;,o 

7375.5.5 
J’Ol’.S 

’4S63.7 
7591 5.8 
”970, q 
410?’. 0 
4  7C  7  7 , 1 
3.3 

’061.1 
4i:’.2 

5153.7 
7’3-t.i. 
10 ’55, 5 
123=  6 . 6 

14757.7 

15400.8 

10463,0 

2161  i.o 
?’56^.1 
2-4613,2 
26664.: 

23715.4 

307^5.5 
3’?  1  7, 6 
’4^6  3^7 

36913.6 
739  7  0 .9 

41022.5 
4  3’ 7  3,  1 

0,0 

2051.1 

4102.2 

5157.3 

3204.4 
1J255.5 

12306.6 
147=7.7 
1640  3,8 
134  =  9  ,0 

20511.6 

22562.1 

24613.2 

26664.3 
’3715.4 
3  QT66 .5 

32817.6 

34363.7 

76919.8 

33973,0 

41022.0 

43077.1 


27774. 0 

’7774.0 
’’774,0 
’’774. C 
27774.C 
27774.0 
’7774.  0 
’7774.  0 
27774.0 
?  ?’?4  •  0 
27774.0 
2  7774.,. 
’7774.0 
27774,0 

’7774,0 

27774.0 
’7774.0 
27774.  ' 
27774,  0 
’7774.0 
’7774.0 
27774.0 

29625.5 
2=625.6 

29625.6 
2=625,6 
296’5,6 

29625.6 
2=625.6 
296’  =  . 6 

29625.6 
2=625.6 

29625.6 
2=6’5.= 
’=625.6 

29625.6 

29625.6 
2=625.6 
2=625.6 

29625.6 
’9625.6 
’9625.6 

29625.6 
2C6’5,  5 
31477.’ 

3 1477.’ 

51477.2 
31477.’ 

3  1477.’ 

31477.2 

31477.2 
31477.? 

3 1477 « ’ 

31477.2 

31477.2 
3_ii*77_1_2 

31477.2 

31477.2 

31477.2 
31477. 2 

31477.2 
314  77,  2 

31477.2 

31477.2 

31477.2 
31477.’ 


SIG 

VAR 

?2»'+r> 

.  34  ==  +  0  3 

,  <  a.  0  r  *.  n  ** 

. <= 7££+ Q3 

,  64  f ’  +  0  2 

ill"**? 

.81  61+0  3 

7  *sc*  -  2 

.  °4  5E+  f  ? 

. ’966+03 

i^i" 

.’C-6E+C3 

•  16+6+ 03 

+ 11 "+  np 

•  1 7’=+  =  3 

C2 

.  75  76  +  0? 

.fo;.***  rp 
,  7^o“4-no 

.49=5+0? 

,47'’+nr 

.7?’"+0? 

« <+* "  ‘ 


•  Z7+eE*  0? 
.  7E  :*E  +  r  z 
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Table  2.  Details  of  the  Statistical  Analysis  (Cont) 
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7.  CONCLUSIONS 


In  Section  6  we  pointed  out  that  the  site  chosen  for  characterization  is  one  for 
which  electromagnetic  scattering  data  is  available.  We  have  also  discussed  how  the 
procedures  of  this  report  were  driven  by  the  intention  to  use  the  results  in  rough 
surface  scattering  calculations.  Briefly,  the  complex  dielectric  properties  of  a 
geological  area,  together  with  knowledge  of  the  mean  height,  variance  in  height, 
correlation  length,  and  PDF  for  the  surface  heights  is  sufficient  for  the  calculation 
jf  the  specular  and  diffuse  reflection  of  a  radar  wave  from  the  surface.  The  re¬ 
sults  of  the  calculations  using  the  statistical  properties  will  then  be  compared  to 
12 

the  data  of  McGarty. 

The  comparison  involves  several  aspects.  Computer  programs  are  being 
developed  which  incorporate  a  number  of  existing  models  of  rough  surface  scatter¬ 
ing  cross  sections,'*'®'  into  -a-more  general  model  for  calculating  the  coher¬ 

ent  and  incoherent  power  detected  by  a  monopulse  receiving  antenna.  The  scatter¬ 
ing  geometry  that  is  used  in  the  model  is  shown  in  Figure  2.  The  terrain  statistics 
of  Table  2  will  be  used  as  one  input  to  this  program  and  the  results  compared  with 
the  experimental  values. 

Once  the  characterization  approach  to  terrain  Scattering  analysis  yields  reason¬ 
able  agreement  with  data,  the  scattering  at  additional  sites  can  be  characterized 
and  these  results  used  to  evaluate  site  terrain  effects  in  detection  and  tracking  of 
low  flying  targets. 


Figure  2.  Reflection  of  Radar  Waves  F  rom  Rough  Terrain 
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Appendix  A 

The  Form  of  the  Jacobian  for  Some  Sets 
of  Related  Variables 

In  the  evaluation  of  multiple  integrals  it  is  often  helpful  to  carry  out  a  change 
of  variables.  This  process  requires  the  determination  of  the  Jacobian  expressing 
the  relation  between  two  sets  of  variables. 

In  particular  we  are  concerned  with  the  sets  of  variables  {Wj,  w2,  . . . ,  wN> 
and  (r,  9  g.  ...»  and  the  set  of  appropriate  connective  relations.  If  we 

can  express  these  relations  such  that: 

(1)  w j  is  a  function  of  r,  Wg,  ....  w^  with  r  alone  varying; 

(2)  w2  is  a  function  of  r,  9  w^,  ....  wN  with  0  j  alone  varying; 


(N)  wN  is  a  function  of  r.  By  6  2»  ...»  9  with  9  j  alone  varying; 

then  the  Jacobian  can  be  expressed  as  the  product  of  the  several  partial  differential 
coefficients  for  those  relations 

9  w,  9w„  3  w„  3  w 


3w3 


N 


57 


N-l 


"See  Edwards,  J.  (1954)  A  Treatise  on  the  Integral  Calculus,  Vol.  II,  Chelsea 
Publishing  Company. 
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As  an  illustration  consider: 

<x,  y,  z)  and  (r,  0,  p)  where 

x  =  r  sin  9  cos  <p,  y  =  r  sin  9  sin  <p,  and  z  =  r  cos  9 
which  we  then  write  as 

x  =  .»/  r2  -  y2  -  z2  z  =  r  cos  (?  y  =  r  sin  9  sin  0 


each  successive  relation  containing  an  additional  member  of  the  second  set  of 
variables.  Then 

9x  9z  9  y  2  .  . 

J=S~F  '  37'3p-r  Sln  0  * 

For  the  particular  sets  of  variables  with  which  we  are  concerned  we  can  estab¬ 
lish  a  similar  series  of  relations: 


(1) 

w2 

=  r  cos  0N_1 

(2) 

W3 

=  r  sin  eN_1 

(3) 

W4 

=  rsin®N_i 

(N-l)  WN  =  r  Sin  ®N-1  sin  6N-2  sin  02  COS  01 

(N)  Wj  =  r  Sin  sin  •  *  *  sin  &2  sin  01' 

This  has  the  proper  form  if  we  rewrite  the  relation  for  Wj  in  its  equivalent  form: 

wi=  Vr2‘ w22- W32-**-  -  wN2  . 


Wj  =  r  sin  sin  •••  sin  9 2  sin 


Then  the  Jacobian  for  our  variables  is 

J  =  <-l)N_1  rN  1  sinN  2  #N_j  sin^"3  si°2  ^3  sin  9  2 
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Appendix  B 

Evaluation  of  the  Zeroth  and  Second  Moment 
Integrals  of  a  Multivariate  Exponential 
Probability  Density 

The  purpose  of  this  section  is  to  present  the  mathematical  analysis  required  to 
evaluate  some  particular  multiple  integrals  involving  sinusoidal  forms.  Specific¬ 
ally,  the  integrations  of  interest  are  those  required  for  the  evaluation  of  the  zeroth 
and  second  moment  integrals  for  the  selected  N-variate  exponential  probability 
density. 

The  zeroth  case,  the  simpler  one,  involves  successive  integrals  of  descending 

powers  of  sine  terms.  The  result  is  also  of  value  in  the  more  complicated  second 

moment  evaluation. 

it  nl  2 

Consider  that  J  sin^  0  60  =  2  j  sin^  6  60  and  that  we  have  the  relation: 
o  o 

1 

B(p, q)  =  /  xp_1  (l-x)q“1  dx 
o 

where  B(p,q)  is  the  Beta  function. 

2 

Then  for  x  =  sin  0  and  dx  =  2  sin  0  cos  0  d@ 
jt/2 

B(p.q)  =2  /  (sin  0)2p'2  (cos2  0>q"1  sin  0  cos  0  60 
o 


and 
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/  sin  9  d0  =  B  ,  |)  . 


For  the  zeroth  moment  case  we  have: 


n 


IT 


7?  IT 

lQ(6)=  J  sinN'2  0N_j  d0N1  /  sinN'3  0N_2  d@N_2  *•*  /  sin203  d03  /  sin  ©2  d02 
o  ‘  o  o  o 


So, 


ye) 


~[B  (V'  i 


*  *>]• 


Now 


b(^±1,  |)  =  r(M+J.)r(|)^r{M+2) 


which  leads  tc 


V*>  = 


N  -2 


<tt>  r~  r  (^1)  r  (|>r  (d 

r(|)  r  (^-i)  •**  r(|)  7 


N  -2 


'2' 


When  combined  with  the  remaining  terms  the  final  value  for  the  total  zeroth 
moment  is: 

2  (jr)N/2Cir(N) 

ITT 


I  = 
o 


y.\1x2----  an  c2w  r  (n/2) 


=  1 , 


The  evaluation  of  the  second  moment  is  more  complex.  Since  the  form  is 


l2  ~ 


^XlX2  *’*  XN 


/  ***  / 


“  "(wi2  ,V,. 

\  *  i  2  N  / 


-*>  \  i  2 


-Cjjy.wj2  *-  w/ 

#  *•  J 


onl/2 


dW  j  •  •  •  d\V*s 
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we  must  make  use  of  both  the  Jacobian  and  the  connective  relations  between  the  two 
sets  of  variables  as  described  in  Appendix  A.  The  new  form  is; 


I2  = 


r  N+l  ~C2r  . 

J  re  dr 


L^A1  A2  ~  AN  ° 

[*  V 

X~  /  sinN  2  0N-1  COS  ®N-1  d®N-l  f  ““  0N-2d0N-2**‘ 

2  o  o 


n  n  2it 

/  sin2  fl3  dS 3  /  sin  02  02  d 0g  / d0x 
o  o  o 


+  XT  /sinN  flN-l  dflN-l  J"  sinN  3  eN-2  008  0N-2  dflN-2  fsin  *  flN-3  dflN-3 
3  o  o  o 


2* 


/  sin  0  2  0  2  d®  2  S  d0l 


4  • 


»  »  *  2j 

+  XT  /si”N(,N-ld0N-JsinN"l9N-2dV2'”  /sin3  02  d02/  cos2  0  d0j 
Wo  o  o  o 


v  n  v  2» 

+  1  /sinN  dfl  j  sin*'1  e  d0N_2  /  sin3  02  d02  /  sin2  0  dfl 

1  o  o  o  o  J 
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In  a  similar  fashion  to  the  previous  case  we  make  use  of  the  relations: 


N+l  C2r  ,  (N+l)l 
Jr  e  dr=  , 

o  C2 

f  sinM  6  d 6  4^  « 


and 


J  sin  ft  cos  ft  dft  s  »  ?>)  - 

o 


These  substitutions  lead  to: 

/  2Cj  (N+1H  \ 

12  ■( 


(0) 


where 


i2«n  = 


|)  •••  B(I.  |)B(4.  |) 


b  (— 4,b(t^'  |>  •**  B»'  !>b4*  4> 

+ — — - j. - - — ”"-r— 

A3 

B  (^i,  |)  B  (|,  4>  B  (^,  |)  •  *  *  B  «i.  4>  B  <4-  |) 


B(‘^i.  yB(|,  4>B(1,  |»B(4.  |) 

XN-1 

b{¥*  1)b{1*  |>***  B(l*  4>b<2.  !>b<!’  I} 

XN 

B  \)  B  (|,  4>  ’  *  *  B  C§.  |>  B  <2.  |>  B  (§,  4> 


T 


1 


1 9(0  )  is  thus  a  complicated  sum  of  products  of  Beta  functions.  It  is  then  necessary 
to  evaluate  this  term.  The  procedure  involves  showing  that  all  the  individual 
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products  of  Beta  functions  are  equal.  The  last  term  (A.  =  A^)  is  a  special  case  where 
none  of  the  integrands  contain  a  cosine  term.  This  will  be  considered  separately. 

In  the  remaining  terms,  there  are  (N  -  1)  integrands,  one  of  which  has  a  cosine  - 
squared  factor  multiplying  some  power  of  a  sine  function  and  the  others  are  ex- 
pressable  as  powers  of  various  sine  functions.  In  each  successive  integrand  one 
additional  sine  -  squared  factor  is  present  compared  to  the  preceding  product.  This 
suggests  that  by  looking  at  the  general  Kth  and  (K  +  l)St  terms  and  showing  equiva¬ 
lence  we  can  demonstrate  that  all  are  identical. 

The  general  term  is: 

u^Di  =  a:+\  B<^,  !>  TT  b  J[  B(^i,  ±> . 

j=i  3=i+l 

If  we  consider  the  differences  in  form  for  the  K**1  and  (K  +  l)s*  terms  we  find  two 
in  each: 

A  „,N-K  3  \  n/N-K-1  14 

AK=  B(— '  2)B{— ’  V 


„,N-K+2  1  \  0/N-K-l  3, 

AK+1  =  B(  2  *  2  B  l  '  2  *  • 


The  remaining  product  terms  are  the^same  in  both  cases.  Thus  the  two  terms 
will  be  equivalent  if  we  can  show  that  A^,  =  AK+1>  To  show  this  we  write  the  Beta 
functions  as  Gamma  functions: 

r  <^>  r  (§)  r  r  <|>  r<|)  r  (i)  r 


r  t5fcK±£)  r  (i)  r  (2L£J)  r  (J, 


r  (I)  T(|>  r$kKz±) 


So  A „  =  A„  .  and,  in  general,  consecutive  terms  will  be  equivalent.  For  the 
special  case  of  the  Nxn  (A.  =  A term,  consider  that  for  K  =  (N  -  1)  the  general 
expression  is: 

R,1  3,  N-2 

i^vi--V  j; 
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while  for  K  =  N 


[i2<0)!n 


(1  i 
2'  2 
T 


N-2 

IT 

j=i  ‘  ‘ 


These  two  expressions  differ  only  in  the  form  of  the  initial  factors.  Since  the 
Beta  function  is  symmetric  with  respect  to  its  arguments  (as  is  clearly  seen  when 
it  is  written  using  Gamma  functions)  these  final  terms  are  also  equivalent  and  we 
have  shown  that  all  the  products  of  Beta  functions  determining  I2(0 )  are  equivalent. 
Since  all  the  Beta  function  products  are  equivalent,  we  then  have: 


N 


*2  = 


'  2  C.  (N+l)l  _  ... 

1  i=l  i 


HE 


N-2 


XN  C2 


■H+Z 


B<: 


b m±.h. 


j=i 


To  simplify  we  expand  the  Beta  function  and  obtain: 

B(^,  I1  IT  j)  »  (tt)N/2  [2r(^)]_1 


(u)N/2  [Nr(N/2)f1 


and  hence: 


N 


2  C.  r<N+2)  £  A .‘l  (v)N/2 
1  i=l  1 

N\IX1X2  *  "*n  C2N+i!  F 


No 


This  completes  the  evaluation  of  the  multiple  integrals  required  to  specify  the 
unknown  coefficients  of  the  multivariate  exponential  probability  distribution  function 
for  use  with  the  data  sets. 
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